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ABSTRACT 

Interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way provide 
the majority of the gamma rays observed by the Fermi Gamma Ray Space Telescope. In addition 
to the gas which is densely concentrated along the Galactic Disk, hydrodynamical simulations and 
observational evidence favor the presence of a halo of hot (T ~ 10 6 K) ionized hydrogen (Hn), 
extending with non-negligible densities out to the virial radius of the Milky Way. We show that 
cosmic ray collisions with this circum-galactic gas should be expected to provide a significant flux of 
gamma rays, on the order of 10% of the observed isotopic gamma ray background at energies above 1 
GeV. In addition, gamma rays originating from the extended Hn halos of other galaxies along a given 
line-of-sight should contribute to this background at a similar level. 

Subject headings: galaxies: evolution - Galaxy: halo - cosmic rays - gamma rays - methods: numerical 



1. INTRODUCTION 



In addition to emission clearly associated with a Galac- 
tic origin, the Fermi Gamma Ray Space Telescope (as 
well as it s predecessors, SAS-2 ; Fichtel et al.||1978{ and 



EGRET; Sreekumar et al. 1998) 



ras identifie d an approx 



imate ly isotropic background of gamma rays ( Abdo et al. 
20101 which is generally considered to be extragalactic 



m nature. For many years, unresolved blazars were con- 



sidered likely to be a dominate contributo r ( Stecker fc 
Salamon] |T996; Kneiske k Mannheim 2005| |A"bazajian et 



al.|201ip , although the lack of observed anisotropy in the 
gamma ray background constrains the blazar contribu- 
tion to be less than about 20% of the observed emis- 



sion (Ackermann et al.||2012[ |Cuoco et al. 2012 



however, Harding & Abazajian 2012" ~. In light of this 



see, 



starburst galaxies are now perhaps the most promising 
class of so urces to generate the ma jority of this observed 



emission (Thompson et al. 2007). Ot her extragalactic 



sourc es, including gamma ray bursts ( [Casanova et al . 
2007[), s hocks associated with large-scale~structure tor 



maHbn (ILoeb fc Waxmanj |2000| |Keshet e t~aL| [2003]) 



ult ra-hig h energy cosmic ray interactions (Kalashev et 



al. 2009|), and dark matter annihilation s or decays (UT 



lio et al. 2002 Belikov & Hooper 20101 have each been 



considered as possible contributors to the observed back- 
ground. 

And while the roughly isotropic nature of this back- 
ground is suggestive of an extragalactic origin, it is diffi- 
cult to exclude the possibility that a significant fraction 
of this e mission is produced more locally. It was sug- 
gested by Keshet et al. (2004), for example, that inverse 
Compton scattering ot high-energy electrons in a halo ex- 
tending out to at least several tens of kiloparsecs (kpc) 
beyond the disk of the Milky Way could potentially ac- 
count for much of the "extragalactic" background. In 
this article, we consider an alternative Galactic contri- 
bution to the isotropic gamma ray background. In par- 
ticular, we consider cosmic ray proton interactions with 
a spatially extended halo of hot and ionized gas. 

This scenario is particularly plausible for a number 



of reasons. Firstly, there is an increasingly compelling 
body of evidence in favor of the presence of an extended 
circum-galactic reservoir of hot (T ~ 10 6 K) and ionized 
gas. Observationally, thi s includes the detection of dif- 
fuse soft X-ray emission flKuntz k Snowden||2000| |Lumb 



et al. 112002 Hickox k Markevitch 200bL the presence ot 



absorption lines in the ultraviolet and in X-ray bands 



consistent with a diffuse background ( Sembach et al 
20031 INicastro et al.||2003l |Rasmussen et a " 



rcg- 



man||2007p, Ha emission and morphological peculiarities 



ot the Magellanic st ream caused by interactions with the 
cir cum-galactic gas (|Weiner k Williams 1996; Putman et| 
aT|2003| |Bland-HawtEorn et al.||2007[ |Besla et al.|2012[), 
th e distorted structure of high velocity clouds ( |Briins et| 



al. 2000 Peek et al 



2007[|Putman et al |[20TT| > , and the 



Hi deficiency of Milky- Way dwarf galaxies as a result of 



ram-pressure str ipping ( Blitz k Robishaw|2000 Grcevich 
k Putman|2009 ) . And while the spatial extent and total 



mass ot this circum-galactic halo are rathe r uncertain 



sizes of ~100 kpc or more are plausible ( e.g., Blitz k Ro 



bishaw||2000l IGrceyich fc Putman|[2009l cf. [Anderson fc' 
Bre gman|| 2010[ ). Simulations and semi-analytic models 



each predict that the circum-galactic halo should extend 
out to the virial radius, and contains a large fraction or 
most of the so-called missing baryons (|Toft et al.||2002[ 
Mailer fc B ullock 2004 Sommer-Larsen|200tit |Rasmussen| 
et al.||2009| |Crain et al.||2010| |Gnedin||20irf 



Secondly, as only a small fraction of the cosmic rays 
produced throughout the history of the Milky Way are 
presently confined to the Galactic disk (see Sec. [3]), one 
expects there to be a spatially extended halo of cosmic 
rays. For reasonable estimates for the diffusion of cosmic 
rays in the outer halo of the Milky Way, we find that the 
rate of interactions between cosmic rays and the circum- 
galactic gas are sufficient to produce a significant fraction 
(on the order of 10%) of the observed isotropic gamma 
ray background. 

The remainder of this article is structured as follows. 
In Sec. [2] we discuss theoretical predictions which fa- 
vor the existence of an extended halo of hot, ionized gas 



2 



around the Milky Way and compare them with observa- 
tional data. In Sec. [3j we discuss the characteristics of 
the cosmic rays in this outer halo, and in Sec. [4] we es- 
timate the gamma ray flux from their interactions with 
the circum-galactic gas. Finally, in Sec. [5] we discuss our 
results and draw conclusions. 

2. THE DISTRIBUTION OF GAS IN THE OUTER HALO OF 
THE MILKY WAY 

While gas in the form of neutral atomic and molec- 
ular hydrogen is strongly concentrated within the disk 
of the Milky Way, a halo of ionized hydrogen (Hn) may 
extend out to the virial radius of the galaxy. In this 
section, we discuss estimates of the distribution of this 
circum-galactic gas, based on a high resolution cosmo- 
logical, hydrodynamical simulatio n run with the ART 
code (Adaptive Refinement Tree; Kravtsov et al.|["l997| 
2002 1 . As the set-up and details ot the physical model- 



mg are similar to that of previous n umerical experiments 



performed by some of the authors (Gnedin & Kravtsov 
2011 Feldmann et al.|2011 ), here we give only a brief re- 



count and highlight the differences between these various 
simulations. 

ART is an Eulerian hydrodynamics + N-body code 
based on the adaptive mesh refinement technique that 
allows one to increase the resolution selectively in a speci- 
fied region of interest, here within five virial radii (i? v ir = 
400 kpc) around a randomly selected M v - lv = 2 x 10 12 M 
dark matter halo in an 8.6 Mpc box. The virial mass and 
radius refer to a region that is 180 times denser than the 
mean densitjF] of the universe. The selected dark mat- 
ter halo is embedded in layers of subsequently lower dark 
matter resolution to reduce the computational cost, while 
still being able t o correctly capture the impact of large 
scale tidal fields ( |Katz||1991J |Bertschingerj|2001[ ) . 

The simulation is started from cosmological i nitial con - 
ditions consistent with WMAP data ( |Spergel et al.|2003[ ): 
h = 0.7, n A = 0.7, Q m = 0.3, n b = 0.043, (implying 
a universal baryon fraction of /j, = f2t/f2 m = 0.143), 
an d as = 0-9- In contrast to the simulations presented 
by Gnedin & Kravtsov (2011 ), this simulation is contin- 
ued self-consistently down to z — 0. 

The simulation includes a photo-chemical network to 
compute abundances of t he various hydrogen (including 
H 2 ) and helium species ( |Gnedin et aLl|2009] [Gnedin fc| 
Kravtsov||2011|). It accounts tor metal enrichment from 



type la and type II supernovae, but not for any additional 
thermal or momentum feedback. Metal-dependent radia- 
tive cooling rates are computed in the optically thin limit. 
Radiative transfer of UV radiation from stellar sources is 



follow ed in the OTVET approximation (Gnedin & Abel 
2001 1. Star formation is based on the amount of molec- 
ular hydrogen present. 

The simulation has a peak spatial resolution of 250 
pc in comoving coordinates. Dark matter particles in 
the high resolution region have a mass of 1.4 x 10 6 M©, 
while the masses of stel lar particles depend on the star 



formation rate, see e.g., Feldmann et al. (2012), but are 
no smaller than 10 3 M Q . 
At z = 0, a massive disk galaxy resides at the center of 

1 Alternatively, M 20 o = 1.7 X 10 12 M© and i? 2 oo = 243 kpc are 
the quantities corresponding to a region that is 200 times denser 
than the critical density of the universe. 
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Fig. 1. — The density of ionized hydrogen (Hn), neutral hydrogen 
(Ht), and molecular hydrogen (H 2 ) as a function of distance to 
the Galactic plane, relative to a location 8 kpc from the Galactic 
center, as found by the simulation described in the text. The inset 
shows the Reynolds layer of Hn above and below the midplane and 
part of the extended diffuse circum-galactic component. Note that 
although the Hi and H 2 densities fall off to negligible levels outside 
of the Galactic plane, a significant density of Hn gas persists to 
the virial radius of the halo. 

the simulated Milky- Way-like halo. It is immersed in a 
low density halo of circum-galactic gas that extends out 
to the virial radius of the halo. In Fig. [T] we show the 
density profile of ionized, atomic and molecular hydro- 
gen as function of the height above the Galactic plane 
centered on a midplane annulus with a radius of 8 kpc 
relative to the Galactic center. Atomic and molecular hy- 
drogen are strongly confined to the Galactic plane and 
drop to negligible densities ~ 1 — 1.5 kpc above the mid- 
plane. A Reynolds layer of ionized hydrogen surrounds 
the neutral components of the interstellar medium and 
extends out to ~ 2 kpc above the midplane. But ionized 
gas is also present at much larger radii. In fact, the sim- 
ulation predicts that ionized circum-galactic gas with a 
density of ~ 10~ 4 cm~ 4 is present out to ~ 100 kpc and 
beyond. 

In Fig.[2j we compare the predictions of this simulation 
to observations. Note that we switch from cylindrical 
profiles at the solar circle in Fig. [I] to spherical profiles 
around the galactic center. The predicted density profile 
of the circum-galactic Hn is consistent with density esti- 
mates and upper limits based on observations. We note 
that the uncertainties in the current observational esti- 
mates are substantial (if constrained at all) , especially at 
large radii (> 100 kpc). Nonetheless, the agreement is 
encouraging. 

Fig. [2] also shows that the shape of the density profile 
of the circum-gal actic gas is different from the one of the 
total matter (cf. Anderson & Bregman 2010). While a 
rescaling of the total matter density profile with a fac- 
tor of ~ 0.05 /b allows to match the observed density 
at ~ 50 kpc, it leads to a severe underestimation of the 
density at larger radii and the overall mass in the circum- 
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of a diffuse X-ray halo around the gian t spiral galaxy 



2011| ). The observa- 
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Fig. 2. — Comparison between our simulation and observations. 
Shown are the total matter density profile as measured in the sim- 
ulation (blue solid line), the total matter density profile of a dark 
matter halo with NFW profile and M 2 oo = 2 X 10 12 M Q and 
R200 = 250 kpc (blue dot-dashed line), the profile of Hn as mea- 
sured in the simulation (red solid line), a fit of a beta-profile to the 
Hn density profile over the range 30-400 kpc (red dotted line), and 
the total matter density profile rescaled by a factor 0.05 /j, (red 
dashed line). Note that, in contrast to Fig. [I] density profiles are 
measured in spherical bins around the center of the galaxy. The 
data points indicate observational estimates of the circu m-galactic 
gas (filled symbols) or upper limits (empty symbols) jWeiner fc| 
Williams 1996 Blitz & Robishaw 2000 Grccvich k, Putman 20Uy 
Quilis & Moore 2001 Stanimirovic ct al. 2002, Anderson & Brcg- 
man||2U10p . The total matter density distribution is well approx- 
imated by an NFW profile, except in the central region of the 
halo (< 20 kpc), where cooling and adiabatic contraction lead to a 
density enhancement. The Hn density profile that our simulation 
predicts is in good agreement with observations. A rescaling of the 
total matter density profile by a factor 0.05 /b matches the Hn den- 
sity at r ~ 50 kpc, but severely underestimates the gas density at 
larger radii and the overall mass in the circum-galactic component. 
Our simulation predicts that the baryons in the circum-galactic Hn 
account for ~ 0.25 /b of the total mass in the halo. 

galactic Hn component. In fact, our simulation predicts 
that the Hn density dist ribution is reasonabl y well de- 



scribed by a beta-profile (Forman et al. 
of more gradually with radius than an i> 



1985 ) that falls 
"FWprofile. Af- 



ter subtracting the Hn mass within the central 25 kpc, 
the predicted mass of the circum-galactic ionized gas is 
2 x 10 9 M Q within 50 kpc, 1 x 10 10 M within 100 kpc, 
3.5 x 10 10 Mq within 200 kpc, and 7.3 x 10 10 M within 
i? vir = 400 kpc. The latter mass corresponds to ^0.25 
/h of the virial mass of the halo and is comparable to 
the combined mass of the observed stellar and neutral 
gas components in the Milky Way. Consequently, the 
baryonic fraction of the Milky Way halo is ~ 50% of the 
universal value /b, i.e., half of the baryons must have 
been expelled from the halo or beenprevented from be- 
ing accreted onto it in the first plac^J 
We note that this result is in line with the detection 

2 If the baryon fraction of the Milky Way halo were close to 
the universal value /b, the circum-galactic gas densities and our 
predictions for the gamma-ray fluxes would need to be increased 
by about a factor 3. 



NGC 1961 (Anderson & Bregman 

tionally derived mass ot the ionized gas within ^500 kpc 
is ^(1 — 3) x 10 11 M Q . This correspond^ to an Hn mass 
fraction of ^0.1 — 0.35 f\>, similar to our prediction of 
^0.25 /b- Interestingly, the picture may be different for 
elliptical galaxies. Recent observations of X-ray halos 
around elliptical galaxies of relatively low halo mass (a 
few 10 12 Mq) indicate that the fraction of baryons in 
their ionized gas halos is significantly higher ~ 0.5 — 0.7 
/b, which together with their large stellar mass fraction 
is sufficie nt to account for essent i ally a ll of the "missing 



baryons" (Humphrey et al. 2011 2012). 



3. AN EXTENDED HALO OF COSMIC RAYS 

In galaxies such as the Milky Way, cosmic rays are 
thought to originate from sources distributed through- 
out the disk, such as supernova remnants. These cosmic 
rays proceed to propagate diffusively through the Galac- 
tic magnetic field, resulting in the formation of a steady- 
state distribution which is present throughout the vol- 
ume of the disk. Measurements of unstable species in the 
cosmic ray spectrum can be used to deduce the timescale 
of how long cosmic rays remain in the disk before escap- 
ing. The 10 Be-to- 9 Be ratio is particularly useful in this 
regard, as 10 Be is the longest lived an d best measured un - 
stable secondary cosmic ray species ( Strong et al.[ 2007). 
These measurements lead to an estimate of the cosmic 
ray escape time of t csc sa 2 x 10 7 years, corresponding 
to only ~0.2% of the star formation history of the Milky 
Way. Therefore only a very small fraction of the cosmic 
rays produced throughout the history of the Milky Way 
are presently confined to its disk. 

Under simple 1-D diffusion (such as the case of diffu- 
sion away from and perpendicular to the disk) , a cosmic 

ray will traverse a distance D ~ \ftK in a time t, where 
K is the diffusion coefficient. Measurements of the lo- 
cal cosmic ray boron-to-carbon ratio imply a value of 
K ~ 10 29 cm 2 /s for ~GeV-TeV protons within the dif- 
fusion zone of the Galactic disk. For such a diffusion 
coefficient, cosmic rays injected since the Milky Way's 
peak in star formation (a few Gyr ago) will be concen- 
trated within ~30 kpc of the disk. Of course, the mag- 
netic fields are expected to be considerably weaker in the 
outer halo than within the disk, potentially facilitating 
much more efficient diffusion. For reasonable estimates 
of the magnetic fields, however, one still expects a large 
fraction of the integrated cosmic ray luminosity of the 
Milky Way to be contained within an extended cosmic 
ray halo inside of the Galaxy's virial radius. A lower 
limit on the spatial extent of the Milky Way's cosmic ray 
halo can be estimated by simply considering the galaxy's 
power output in cosmic rays. Presently, the total energy 
in cosmic rays throughout the volume of the disk is ~ 
1 eV/cm 3 x tt(200 pc)(15 kpc) 2 ~ 4 x 10 66 eV. In com- 
parison, assuming a constant cosmic ray injection rate, 
the total energy injected into cosmic rays throughout the 
history of the Milky Way is ~ 5 x 10 40 erg/s x 10 10 yr 

3 The inclination and hence total mass of NGC 1961 are rather 
uncertain. We use the updated inclination estimate i ~ 65° 
| |Combes et al. 20091, which implie s a circular velocity of ~ 300 
km/s at 34 kpc jrlaan et al.||2008| and thus a halo that is 2-2.5 
times more massive tnan tne halo of the Milky Way. 
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- 1.6 x 10 58 erg - 10™ eV, which is about 2500 times 
higher than is currently confined to the disk. Alterna- 
tively, one can consider the ratio of the duration of the 
Milky Way's star formation history to the escape time of 
cosmic rays in the disk, (10 10 yr)/(2 x 10 7 years)~ 500. 
According to either estimate, the vast majority of cosmic 
rays generated by our Galaxy are no longer contained it 
its disk. If we also take into account variations in the 
star formation rate (which was a factor of a few higher 
than present rates between roughly 2 and 6 billion years 
ago), we find that only about ~0.08% of the cosmic rays 
created over the history of the Milky Way are currently 
present within the disk. If we (unrealistically) assume 
that the cosmic rays are distributed with a uniform den- 
sity equal to the value in the disk, the population of cos- 
mic rays would fill a sphere approximately ^35 kpc in 
radius. As measurements of secondary-to-primary ratios 
tell us that the cosmic ray density must drop off signifi- 
cantly beyond a few kpc from the disk, we conclude that 
this halo of cosmic rays must be at least ^50 kpc, and 
possibly much larger, in radial extent. 

If the magnetic fields in the outer halo are on the order 
of ~40 2 times smaller than those in the disk (or that the 
diffusion coefficient is 10 2 times larger than in the disk), 
for example, the majority of the cosmic rays produced 
throughout the history of the Milky Way will remain 
within a few hundred kpc of the disk, corresponding to 
cosmic ray densities in the outer halo on the order of 
0.001-0.1 of that found in the disk. Through interactions 
with the extended halo of ionized gas described in the 
previous section, these cosmic rays can provide a flux of 
gamma rays which constitutes a significant fraction of 
the observed isotropic gamma ray background. 

4. THE CONTRIBUTION TO THE ISOTROPIC GAMMA 
RAY BACKGROUND 

To calculate the gamma ray flux from cosmic ray in- 
teractions with gas in the outer halo of the Milky Way, 
we adopt a simple diffusion model: 



° dN f(x,t,E p ) = V-[K(E p )V^(x,t,E p )]+Q(x,t,E p ), 



dt dE, D 



where K is the diffusion coefficient, and dN p /dE p de- 
scribes the distribution and spectrum of cosmic rays. Q 
is the source term, which describes the injection rate and 
spectrum of cosmic rays from the disk. As the escape 
time from the disk is small compared to the relevant 
timescales of the problem, we simply assume that cos- 
mic rays immediately escape the disk and diffuse out- 
ward according to the K value chosen for the outer halo. 
We adopt an energy dependence of the source term given 
by Q(E p ) cc E~ 2A (consistent with the observed cosmic 
ray spectrum, after accounting for diffusion), and have 
normalized the source term to produce the observed local 
density of cosmic rays. We also adopt a time dependence 
in the source term intended to reflect variations in the 
star formation rate of the Milky Way: 



Q(o) 



l + t/(lGyr) 
3 

v 3-0.5(t-6Gyr) 



if t < 2 Gyr, 

if 2 Gyr < t < 6 Gyr, 

if 6 Gyr < t < 10 Gyr. 
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Fig. 3. — The contribution to the high latitude gamma ray 
backgr ound (as measured by the Fermi Gamma-Ray Space Tele- 
scope; Abdo ct al. 2010) from cosmic ray interactions with ion- 
ized hydrogen (-tin) m the outer halo of the Milky Way. Here, 
we have adopted the gas density shown in Fig. [I] and propagated 
cosmic rays with a diffusion coefficient in the outer halo given 
by K(E P ) = 1.2 X 10 29 cm 2 /s(E p /GeV) - 33 (top) and K(E P ) = 
4 X 10 30 cm 2 /s (Ep/GeV) 33 (bottom). We expect a contribution 
at a similar level from gamma rays that are produced by the same 
mechanism in the Hjj halos of all galaxies along a given line-of- 
sight. 

We consider a diffusion coefficient that is constant 
throughout the outer halo and with an energy depen- 
dance given by K(E P ) = K Q E° 33 . We find that for a 

value of Kq = 1.2 x 10 29 cm 2 /s, the cosmic ray halo 
extends out to ~60 kpc, and the cosmic ray density in 
the region surrounding the disk is about 30% of the lo- 
cal cosmic ray density. Much smaller values of K are 
likely inconsistent with measurements of local primary- 
to-secondary cosmic ray species. We also consider a sub- 
stantially larger diffusion coefficient of Kq = 4 x 10 30 
cm 2 /s, which leads to a cosmic ray halo that extends out 
to several hundred kpc. 

From the cosmic ray distribution and the distribution 
of ionized Hh gas as shown in Fig. [T] we calculate the 
spectrum of gamma rays produced per volume from pion 
production: 



dN 7 
dK 



= 2 



W dE„ JeT- 



(2) 



where dN^/dE^ is the spectrum of neutral pions pro- 
duced in cosmic ray-gas collisions: 



dN„ 
dK 



= 47T riH 



dE p J p (E p )-—^(E n ,E p ). (3) 



(1) 



Here, «h is the number density of gas (as shown in Fig.fT]) 
and Jp(Ep) is the cosmic ray intensity (per energy). For 
useful parameterization s of the differential cr oss section 
for pion production, see Blattnig et al. (2000). 

In Fig. [31 we show the contribution to the gamma ray 
background at high Galactic latitudes from cosmic ray 
interactions with the extended halo of ionized gas. Re- 
sults are shown for two choices of the diffusion coefficient, 
which determine how far the cosmic rays have propagated 
after escaping the disk. For this range of diffusion coef- 
ficients, we find that between 3 and 10% of the isotropic 
emission observed by Fermi originates from these inter- 
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actions. 

The spatial extent of the cosmic ray halo represents the 
most significant uncertainty in our calculation. The two 
curves shown in Fig.|3jrepresent current cosmic ray distri- 
butions which fall tohalf of their density (not including 
the density within the disk itself) by 60 kpc and by 360 
kpc, respectively. For these two cases, half of the ob- 
served gamma rays originate from within 18 and 38 kpc, 
respectively, leading to an approximately isotropic an- 
gular distribution of gamma rays, not dissimilar to that 
predicted in t h e inve rse Compton scenario described by 
|Keshet et al.| ( |2004[) . We consider these two cases to 
represent a reasonable range of possibilities, although 
gamma ray fluxes higher or lower by a factor of a few 
are not implausible. 

The overall luminosity of gamma rays originating from 
cosmic ray interactions with circum-galactic gas is po- 
tentially substantial. The large number of cosmic rays 
in the Galactic halo (see Sec. [3| compensates for the 
low density of the circum-galactic gas. We can obtain 
an order of magnitude estimate by integrating the cos- 
mic ray injection rate over the last 10 Gyr and assuming 
that most of the cosmic rays remain within the virial ra- 
dius of the Galaxy, where they interact with the circum- 
galactic gas. Under these assumptions, we obtain a 
gamma-ray luminosity of L 7 (> 100 MeV) ~ 9 x 10 38 
erg s -1 x ((nH)/10~ 4 cm~ 3 ), where (tih) is the typical 
circum-galactic gas density encountered by cosmic rays. 
For (rifj) ~ 10~ 4 cm~ 3 , the gamma-ray luminosity from 
the circum-galactic gas is comparable with that induced 
by cosmic rays inter acting with the inte rstellar medium 
in the Galactic dis k (Strong et al. 2010) or that of M31 
( [Abdo et al.1[2010 ). 



The gamma ray luminosity from such interactions in 
the extended gaseous halos of other galaxies are also ex- 
pected to contribute to the observed isotropic gamma ray 
background. We can use the estimate for the circum- 
galactic gamma ray luminosity given in the previous 
paragraph to judge whether this contribution is likely 
to be significant. We scale the gamma ray luminosity of 
a given galaxy with its stellar mass, which is an indicator 
of the total number of cosmic rays injected into the halo 
of ionized gas. Then, using th e redshift-depende nt stellar 
mass to halo mass rela tion by Yang et al. (2011 ) and the 
halo mass function by Klypin et al. ( |2011 1 , we estimate 
the contribution to the local gamma ray flux by integrat- 
ing over all redshifts and halo masses. We note that only 
galaxies below z ~ 1 and with M vir ~ 10 11 - 10 13 M Q 
are particularly relevant for this calculation. Galaxies in 
lower mass halos do not contain a sufficient number of 
cosmic rays to contribute significantly, and galaxies that 
reside in very massive halos are too rare. Overall, we 
find that the total contribution from extended gaseous 
halos of all galaxies should account for a few percent of 
the isotropic gamma ray background. This contribution 



is thus similar to that originating in the ionized halo of 
our own galaxy. 



5. DISCUSSION AND CONCLUSIONS 

The Fermi Gamma Ray Space Telescope has revealed 
the presence of very li ttle anisotropy in the diffuse 



gamma ray background (Ackermann et al. 2012 Cuoco 



et al. 2012), leading to the disfavoring of blazars and 



other rare but luminous gamma ray sources as the dom- 
inant contributors to this background. Instead, it seems 
that more faint and numerous sources, such as starburst 
galaxies, are more likely to make up the bulk of the ob- 
served emission. 

In this paper, we have considered an alternative contri- 
bution to the observed isotropic gamma ray background. 
In particular, we have shown that the presence of an ex- 
tended halo of ionized gas (Hn) surrounding the Milky 
Way, as favored by both simulations and observations, 
can provide an important target for cosmic rays that 
are no longer confined to the disk, leading to a signif- 
icant contribution to the high-latitude, diffuse gamma 
ray flux. Using a distribution of Hn gas derived using a 
high resolution hydrodynamical simulation and reason- 
able estimates for the distribution of cosmic rays in the 
outer halo of the Milky Way, we estimate that pion pro- 
ducing interactions taking place within the surrounding 
several tens of kiloparsecs from the Galactic disk may ac- 
count for on the order of 3-10% of the isotropic gamma 
ray background as observed by the Fermi Gamma Ray 
Space Telescope. In addition, the corresponding gamma 
ray emission originating in the Hn halos of galaxies along 
a given line-of-sight should contribute at a similar level. 

In addition to the gamma ray background under con- 
sideration here, it is interesting to contemplate whether 
an extended halo of circum-galactic gas may be con- 
nected in any way to the significan t excess of isotropic 



served by several groups (Haslam et al.||1981 


Reich & 


Reich 


1986 


Roger et al.|1999; Fixsen ct al. 2011 


Guzman 


et al. 


2011 


I. In particular, one could imagine energetic 



collisions with the circum-galactic gas, resulting in the 
production of an approximately isotropic background of 
radio synchrotron. With little information pertaining to 
the magnetic fields present in the region of the outer halo, 
however, it is difficult to make quantitative estimates of 
such a signal. 
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